FERMILAB-PUB-01 /035-T 
HUPD-0104 
RCNP-Th01005 
UT-CCP-P-102 
hep-lat/010302^ 



0(a)-improved quark action on anisotropic lattices and 
perturbative renormalization of heavy-light currents 

Junpei Harada", Andreas S. Kronfeld^''^, Hideo Matsufuru'', 
Noriaki Nakajima'' and Tetsuya Onogi'^ 

Department of Physics, Hiroshima University, Higashi-Hiroshima 739-8526, Japan 
^ Theoretical Physics Department, Fermi National Accelerator Laboratory, 

Batavia, Illinois 60510, U.S.A. 
^Center for Computational Physics, University of Tsukuba, 
Tsukuba 305-8577, Japan 
'^Research Center for Nuclear Physics, Osaka University, Ibaraki 567-0047, Japan 

March 22, 2001 
Abstract 

We investigate the Symanzik improvement of the Wilson quark action on anisotropic 
lattices. Taking first a general action with nearest-neighbor and clover interactions, we 
study the mass dependence of the ratio of the hopping parameters, the clover coefficients, 
and an improvement coefficient for heavy-light vector and axial vector currents. We show 
how tree-level improvement can be achieved. For a particular choice of the spatial Wilson 
coupling, the results simplify, and 0{mQar) improvement is possible. (Here toq is the 
bare quark mass and a-r the temporal lattice spacing.) With this choice we calculate the 
renormalization factors of heavy-light bilinear operators at one-loop order of perturbation 
theory employing the standard plaquette gauge action. 



1 Introduction 

The anisotropic lattice has become an important tool in lattice QCD simulations. With a small 
temporal lattice spacing ar one can more easily follow the time evolution of correlators, while 
keeping the spatial lattice spacing Oo- comparatively modest [ |l|. This approach is especially 
effective when the signal-to-noise ratio deteriorates quickly, as, for example, in the case of 
glueballs [ ^ . The better signal-to-noise ratio is beneficial also for heavy quark systems [ ^ . In 
addition, it is hoped that the anisotropy can be exploited to reduce lattice artifacts [Q, which 
are a special concern with heavy quarks. 

In current work on heavy quarks, lattice artifacts are controlled with non-relativistic QCD 
(NRQCD) and heavy-quark effective theory (HQET). This is done either a priori, by discretiz- 
ing the NRQCD action [ or a posteriori, by using the effective theories to describe lattice 
gauge theory with Wilson fermions [ ^, ^. These strategies are possible because the typical 
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spatial momenta in heavy quark systems are much smaller than the heavy quark mass. Heavy 
quarkonia have momenta p ~ mQV and k ~ niQv'^, where v ~ 0.1-0.3 is the heavy-quark 
velocity; heavy- light hadrons have momenta only of order Aqcd- In these approaches one is 
left with discretization effects of order (Aqcd a)" from the light quarks and gluons and of order 
a\{p/mQ)'^ from the heavy quark. 

The method of Ref. [ |6| smoothly connects to the usual continuum limit, so one can, in 
principle, reduce discretization effects to scale as a power of the lattice spacing a, but only 
by making a too small to be practical. Klassen proposed using anisotropic lattices with the 
anisotropy ^ = a^ja^ chosen so that raqa^ and pa„ are both small [ ^. Clearly, this proposal 
works only if p is smaller than mg, as in the approaches based explicitly on heavy-quark 
theory. It also works only if renormalization constants have a smooth limit as moflr ~^ 0, 
where ttiq is the bare quark mass. In particular, one would like to be able to expand the 
renormalization constants in powers of moflr even when m^ar, ~ 1. Then it may be possible to 
adjust the improvement parameters of the lattice action (and currents) in a non-perturbative, 
mass-independent scheme [ ^, ^. If, on the other hand, m^aa dependence appears in an 
essential way, then one would be forced back to a non-relativistic interpretation, as explained 
for isotropic lattices in Refs. [ |6|, |^. 

To our knowledge there is no proof that cutoff effects always appear as powers of m^ar- In 
this paper we try to gain some experience by calculating the full mass dependence of several 
(re) normalization constants, first at tree level and then at one-loop in perturbation theory. We 
focus on the Fermilab action [ which is the most general action without doubler states, 
having different nearest-neighbor and clover couplings in the temporal and spatial directions. 
This action has been applied on anisotropic lattices to the charmonium system [ ^, |l^, 11, 12| , 
as have some actions with next-to-nearest neighbor interactions. The self-energy has been 



calculated at the one- loop level in perturbation theory [ 13|. 

In the numerical work on charmonium, two different choices for tuning the spatial Wilson 
term have been made. One choice is that of Refs. [ where = 1/,^. Another choice is 

that of Refs. [ 11, 12|, where = 1. In the first part of this paper, we study improvement 
conditions for these two choices, as a function of the heavy quark mass. (In a perturbative 
calculation more generally improved actions with rj = ^C^Vs also have been considered [ 13 1.) 



By studying the full functional dependence on ^ and moar, we can test whether moao- appears 
in an essential way. We find that the limit of small rriQar is benign at the tree-level only for 
the first choice, = 1/^. For the other choice, = 1, the continuum limit is reached only for 
moao- <C 1. 

It turns out that with the first choice (r^ = 1/^) two of the improvement parameters vanish 
at the tree-level as rriQar — > 0. This simplifies the one-loop analysis, so in the second part of 
the paper we concentrate on this choice. This calculation has two purposes. The first is to 
study cutoff effects of the renormalization coefficients and to test at the one-loop level whether 
they still appear only as powers of rriQar ■ The second is for phenomenological applications to 
heavy-light matrix elements. Even if a non-relativistic interpretation is necessary, anisotropic 
lattices are a good method for reducing the signal-to-noise ratio [ |^ . 

This paper is organized as follows. Sec. § describes the quark field action and discusses its 
parametrization in detail. In Sec. ^, the expression for the one-loop perturbative calculation is 
given. The numerical result for these perturbative constants are presented in Sec. The last 
section is devoted to our conclusions. We give the Feynman rules in Appendix ^ and explicit 
expressions for the one-loop diagrams in Appendix 
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2 Anisotropic quark action 



This section describes the actions with Wilson fermions [ 14| on anisotropic lattices. We denote 
the renormalized anisotropy with ^, and the spatial and temporal lattice spacings with a^- and Ur 
respectively: 

Oct = COr- (2.1) 

These lattice spacings would be defined through the gauge field with quantities such as the 
Wilson loops or the static quark potential. We therefore consider ^ to be independent of the 
quark mass. 

2.1 Quark field action 

Following Ref. let us introduce an action with two hopping parameters [ 14 1 and two clover [ 



15[ coefficients, 

n 

-Ks^[{rs - 7i)Un,i1pn+i + {rg + 7i)t^I_i,iV'n-i] (2.2) 
i 

+ ^rsCBKs ^ eijk(yijBn,k'4^n + iCRl^s X! ^4iEn,i1pn ■ 
i,j,k i 

This is the most general nearest-neighbor clover action. Note that the notation is slightly 
different than in Ref. [ ^; cb of Ref. [ ^ corresponds to rgCB in (|2.2D. 

It is helpful to change to a notation with a quark mass. We rescale field by 

3/2 

^Pr, = ^V(x), (2.3) 

and introduce the bare mass 

moar = - [rt + 3r,C], (2.4) 

with = Ks/nf Then one can rewrite the action as 

S = o,a3^V^(x) [mo + i(n + 74)1)4 -i(rt- 74)1)+ + ia,eV,CA(3) 

X 

-^arfrsCBC^-B - ^ar^CECa-E] ^{x). (2.5) 

The covariant difference operators D^, D and A^^), and the fields B and E are defined as in 
Ref. [ ^ , except that the lattice spacing a is replaced by Ot or in the obvious way. 

The action S has six parameters mo, r^, C, cb, and ce- Two are redundant and can be 
chosen to solve the doubling problem [|^. In particular, we choose = 1 to eliminate doubler 
states. We then rename = r, but discuss how to adjust it below. The other four parameters 
are dictated by physics. The bare mass is adjusted to give the desired physical quark mass, 
and Cb, and ce are chosen to improve the action. 

Following Ref. [ |1 we also consider a rotated field 

^f(x) = e*^i"-/2 ^ a^dij-D] V'(x), (2.6) 
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where Mi is the rest mass, defined and given below, and di is an improvement parameter. This 
field is convenient for constructing heavy-light bilinears 



K"" = ^Tm^', (2.7) 



< = r7,l5^\ (2.8) 

which, at the tree-level, are correctly normalized currents for all rriQar- Beyond the tree-level 
one may add dimension- four terms to these currents, and one must multiply with suitable 
renormalization factors. 

The renormalization factors and the improvement parameters C,, cb, ce, and di must, in 
general, be chosen to be functions of moar and the anisotropy ^. Below we shall give the full 
mass dependence to check whether, for small moUr, power series such as 

C(e, moar) = at 0) + moarC'it 0) (2.9) 

can be admitted. If moa^ = (,moar enters into the full mass-dependent expression, this series 
would not be accurate when moQa ~ 1. In the past [ ^ the behavior in ( p.Sj ) was implicitly 
assumed. If the expansions of the form ( |2.9| ) do work, then for full 0{a) improvement one must 
adjust C(^,0), C'itO), cb(^,0), and Cij(^,0) in the action, and Zjj,(^, 0), Z'j^{^,0), and di(C,0) 
of the currents Jr = V"^, A^. 

From ( |2.5| ) one can see that conditions for the improvement coefficients can be obtained by 
simply replacing 

c ^ ec, (2.10) 

Ts ^ Cr, (2.11) 
CB ^ ^rcB, (2.12) 
CE CE, (2.13) 

in formulae in Ref. [ ^ . For example, the energy of a quark with momentum p is given by 

where pi = 2a~^ sm{^piarr) and Si = sin(pjao-). For small momentum E"^ = Ml -\- 
p^Mi/M2 + 0{p^a'^), where the rest mass Mi and kinetic mass M2 are 

MiQr = In(l-hmoa^), (2.15) 

1 - e( — S — T+^^V (2.16) 



M2ar \moar{2 + moar) 1 + modr 

To obtain a relativistic quark one sets the rest mass and kinetic mass equal to each other. This 
yields the condition 



/ £,rmoar(2 + moar) \ moa^(2 mpa^) grmoa^(2 mpg^) 
~ V V 4:{l + moar) ) 21n(l-hmoa^) 4(1 moar) ' 

which can be read off from Ref. [ ^. Matching of on-shell three-point functions yields the 
conditions 

CB = 1, (2.18) 

(CC)^-l , , gVmoar(2 + moar) 
Ce = -, tH 1 -, 775 [z.Lu) 

moar{2 + moar) 1 + moar 4(1 + moar) 
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on the clover coefficients, and 



mQar{2 + mofflr) 
eC[2(l+moa,)2 



2M2ar 

(,rmoar{2 + moa,-) 



2eC(l + moaO] 



(2.20) 



(2.21) 



2moar(2 + moaT-)(l + moa^) 

on the rotation parameter. These tree- level formulae ( 2.14D -( 2.21[ ) have been obtained inde- 
pendently by M. Okamoto [ 16|. We see that essential dependence on moa„ = S,mQar indeed 
may arise, depending on how r is tuned. 

From the energy ( p. 14 ) one can also find the energy of states at the edge of the Brillouin 
zone. The energy of a state with n components of p equal to ir/aa is 

Enttr = ln(l + moar + 2nrC). (2.22) 

Although there is some freedom to choose r, discussed below, one still wants to keep En and 
Ml well separated. 

For small niQar the interesting Taylor expansions are 

C = C-^ + ^{r^-r)moar + Oi{moarf), (2.23) 

CE = ^{l+^r) + 0{moar), (2.24) 

di = I{l-Cr) + O{moar). (2.25) 

With the mass dependent factor in ( p.6| ) there is no mass dependence at the tree-level in the 
currents' normalization factors. 

Let us now discuss the choice of the redundant coefficient of the spatial Wilson term r. 
Two choices have been used in numerical calculations: 

(i) r = 1/^ [ ^] . This is a natural choice because then the mass form of the action takes 
a symmetric-looking form, without In the small moflr limit, the tree-level improvement 
parameters become 



c(e,o) 
c'(e,o) 

di(e,o) 



-1 



(2.26) 
(2.27) 
(2.28) 
(2.29) 
(2.30) 



A key advantage is that rriQa^ does not appear; the continuum limit is reached for 
small mQa-r- Furthermore, both C'('^)O) ^-iid (ii(^,0) vanish at the tree-level, which is 
especially helpful in one-loop calculations. A disadvantage is that with r = C, = 1/^ the 
energy splitting between physical states and states at the edge of the Brillouin zone is not 
large. One can circumvent this problem to some extent by choosing appropriate cutoffs 
and anisotropy [ 0] . 

(ii) r = 1 [ 1^, 0, [l^. Now all hopping terms in (|2.2| ) have projection matrices ^(1 =F 7/^); 
and the anisotropic nature appears only in C, and cb ce- But now, if one considers 
what happens to the conditions when moUa ^ 1 while niQar <C 1, then 



CEitmoar) 
di{C,moar) 



1 



1 + C (l + ^moaa + limoao 



16moar 



(2.31) 
(2.32) 

(2.33) 
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keeping terms of order (moOo-)^ and moaT^ioOo-- Clearly, the continuum limit sets in only 
when rriQaa ^ 1- Even then C and di are non-zero already at the tree-level. An advantage 
is that the splitting between the physical states and the edge of the Brillouin zone is larger 
than in case (i). 

In passing we mention that Refs. [ ^, 11, ^] take 



^ _ 1 2 moar 



(,2 + ^moflr 



(2.34) 



which agrees with the Taylor expansion ( 2.23 ), but does not agree with the full mass 
dependence (2.17). The denominator of this expression also is of the form that reaches 
the continuum limit only when rriQa^ <C 1. 

It is instructive to examine the difference between the two conditions on by considering 
the full mass dependence of (. Figure Q plots the right-hand side of ( 2.17| ) against moo- := moaa-, 
for several values of ^ and the two choices r = 1/^ and r = 1. The mass in spatial lattice units, 
motj is chosen not because it is a natural variable, but because one usually would first fix the 
spatial lattice spacing Oo- so that pua is small enough, while moflo- ~ 1. One would then choose 
the anisotropy to make mQUr small. For example, let us consider the charmed quark on a 
lattice with = 1 GeV. The quark mass in spatial lattice units is moo- — 1.2, so if = 4, 
then ruQar — 0.3, which seems small. For r = 1/^ one finds C,{^,rnQar) = C(4, 0.3) ~ 0.26, 
which is only 4 percent larger than C(4, 0) = 0.25. In this sense, moa^ — 0.3 is small. On the 
other hand, for the choice r = 1, C(^, moflr) = C(4, 0.3) ~ 0.20, which is 20 percent smaller 
than (^(4,0). Even worse, the Taylor expansion ( p. 23 ) estimates only 0.14. 
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Figure 1: The tree- level relation of C, with the quark mass in the spatial lattice unit, moo-, for 
various ^. The solid and the dotted lines are for choices (i) and (ii) respectively. 
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Thus, only with the choice r = 1/^ does it seem possible to approximate the improvement 
coefficients by the small mQa-r limit. With this choice it seems possible to treat charmed 
quarks, without appealing to the heavy-quark expansion, at accessible spatial lattice spacings 
and anisotropy around 3-4. Lattice artifacts appear under control and there is probably enough 
room between Mi and the energies at the edge of the Brillouin zone, En-, to accommodate the 
lowest excitations of the D meson. On the other hand, it seems that reasonable choices of 
0(j and ^ do not exist for treating the h quark: nribar remains big, requiring a non-relativistic 
interpretation along the lines of Refs. [ ^, ^. 

The choice r = 1/^ requires no tree- level rotation for the quark field. This is a great 
simplification for one-loop renormalization. Then the quark and anti-quark field operators are 
multiplied by the factor exp(Miar/2) = 1 + \Miar + ©((Mia,-)^). With the choice r = 1 one 
would have to include the rotation term for a consistent 0{ar) calculation. In the rest of this 
paper, we therefore focus on r = 1/,^. 



3 One-loop Renormalization 

To carry out one-loop perturbative calculations, we must specify the gauge field action as well 
as the quark action. We begin this section with the gauge field action and remark on the gauge 
couplings, which, in general, differ for the spatial and temporal components of the gauge field. 
Feynman rules required at one-loop level are summarized in Appendix 

The self-energy at the one- loop level is represented by the diagrams in Fig. |3|(a)-(b). We 
calculate, as a function of mga,-, the one-loop contribution to the quark rest mass and wave 
function renormalization factor. These quantities require the self-energy and its first derivative 
with respect to the external momentum p4, evaluated on the mass shell (p4,p) = (iMi,0) [ 
17|. By obtaining the full mass dependence, we can check how the one-loop corrections behave 



for rriQacr ~ 1 and niQar small. We also discuss mean field improvement of the self-energy. In 
the past, the full mass dependence of the one-loop quark self-energy has been obtained for the 
Wilson action on isotropic lattices [ for the Fermilab action on isotropic lattices [ 0] , and 
for several improved actions with rt = ^'^C^s on anisotropic lattices [ 13|. 

We also discuss vertex corrections at the one-loop level, shown in Fig. ^(c), and present 
matching factors for the vector and axial vector currents. We again obtain the full mass 
dependence first, and use it to study the practical situation with moaa- ~ 1 and moar small. In 
the past, the full mass dependence of the one-loop quark vertex functions has been obtained on 
isotropic lattices for the Wilson action [ O] , the clover action [ |l^ , and the Fermilab action [ 



20|. 





P p + k P 




p p 

(a) (b) 

Figure 2: Feynman diagrams for the quark self-energy (a) and (b), and for the vertex correc- 
tion (c). 
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3.1 Gauge field action 

We use the standard Wilson gauge action on the anisotropic lattice [ |l] : 



'S'gauge — PO'rO'cr ^ ] 



3 ^ 3 



i<j=i i=i 



(3.1) 



where Ufj^u denotes parallel transport around a plaquette in the ijlv plane. The bare anisotropy 
7G coincides with the renormalized isotropy ^ at the tree-level. In gauge field theory with 
Nc colors, the coupling /3 is related to the usual bare gauge coupling qq hy j3 = 2Nc/gQ. 

There is a subtlety in the gauge coupling, because the temporal and spatial gluons have 
different couplings [ |l| . One can rewrite (5 and -yc as 

/3cr = ~57 77 = ) (3-2) 

= ^|^ = /57G, (3.3) 

where and are couplings for spatial and temporal gluons, respectively. Although at the 
one-loop level g'^ and g^ need not be distinguished, it is convenient to separate the results for 
spatial and temporal parts. To improve perturbative series, it is crucial to use renormalized 



couplings, defined at the momentum scale typical for the process under consideration [ 22 
These couplings, and therefore the scales, could be defined separately for spatial and temporal 
gluons. With this end in mind, we shall present the coefficients of g"^ and g'^ separately. 

3.2 Rest mass renormalization 

The relation between the rest mass to the self-energy is [ 17] 

gMia, = 1 + ruQar - tr[P+S(iMi, 0)]a^, (3.4) 

where P+ = (1 + 74)72 and the self-energy S(j)4,p) is the sum of all one-particle irreducible 
two-point diagrams. The formula (3.4) is valid for all masses and at every order in perturbation 



theory. Since it is obtained from the pole position, the rest mass is infrared finite and gauge 
independent at every order in perturbation theory [ |2l|]. We write the perturbation series as 

00 

S(ip4,p) = ^<7''sW(p4,p;mo), (3.5) 
1=1 

where we explicitly specify the bare quark mass ttt-q. The quark is massless (Mi = 0) when the 
bare mass is tuned to 

moc = tr[P+S(0,0;moc)]. (3.6) 

It is more convenient to introduce a subtracted bare mass Mq = uiq — moc, which vanishes for 
a massless quark. Then the formula for the rest mass becomes 

gMia, = l + Moa^-tr[P+S(iMi,0;Mo)]a^, (3.7) 

where 

i;(p4,p;Mo)] = i;(p4,p;mo) - moc- (3.8) 
In developing the perturbation series, now Mq is treated independently from g^ . 
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The perturbative series for Mi is 



From 



M 



Ml 



1=1 



a^hn{l + Moar), 
tr[P+SW(iMi,0;Mo) 
1 + Moar 



(3.9) 

(3.10) 
(3.11) 



In evaluating T,^^'^ one may disregard the distinction between Mq and mg, because moc starts 
at one-foop order. To show the mass dependence it is convenient [ ^] to introduce the multi- 
pHcative renormahzation factor Zmi defined by, 



MiUr = tanhMj a^. 



(3.12) 



From = ^a^/ tanh Aff ' a-r one can then remove the anomafous dimension by writing 



^[1] 

'Ml 



Cf 



1Q-7T' 



ln{Mfa: 



(3.13) 



Numerical results for M^^^Qt and c^^^ are given in Sec. ^. 



3.3 Wave function renormalization 



The all orders formula for the wave function renormalization factor is [ 17] 



Z. 



-1 



tr[P+S(iMi,0;Mo)]a^ 



where 



1 dYi 

tip4,p;Mo) = -— (p4,P;Mo). 



In view of the mass dependence, we write 



il + Moar)Z2 = l + Y.g'''zS\ 



(3.14) 
(3.15) 

(3.16) 



1=1 



so that the are only mildly mass dependent. This definition of Z2 is slightly different from 
that of Ref. [ 0] for / > 1. 

The wave function renormalization factor is infrared divergent and gauge dependent. There- 
fore we express the one-loop term as 



= (l + MoaO-itr[P+(sW+sW)], 
=: lyW+LW, 



(3.17) 
(3.18) 



where W^^"^ and L'^^ are the infrared finite and singular parts, respectively. The infrared 
divergence does not depend on the ultraviolet regulator, so it is the same as in the continuum 
theory. We define L'^I by the continuum expression. For a massive quark (A^ <C rn^ <C 07^) 
Feynman gauge, 



L 



[1] 



Cf 
167r2 



(3.19) 
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which we use for the heavy quark. In particular, wji^^ is defined by combining ( p.l8| ) and ( 3. 19] ) 



with m = . For a massless quark (m^ = 0, <C Or^)) the mass singularity seen in ( 3.19| ) 
can still be regulated by the gluon mass. In Feynman gauge, 

LP = ^HalX% (3.20) 

which we use for the light quark. Here Cp = — ^/2Nc [= 4/3 for SU(3)]. Because the 
infrared and mass singularities have been subtracted consistently, we should (and do) find 
limmo^o wj^^ = wl^K Numerical results for wj^'^ and VF/^' are in Sec. ^. 

3.4 Mean field improvement 



Mean field improvement [ 22 1 has been employed extensively in Monte Carlo work to improve 
tree-level estimates of couplings. The approximation works, because most of the one-loop 
coefficients can be traced, via tadpole diagrams, to a mean field. On the anisotropic lattice, 
the mean field values of the link variables can be defined individually for the temporal and 
the spatial links. We denote them by Ur and respectively. Then mean field improvement is 
achieved by replacing the link variables with [ ^, 10, |ll], 

Ui^U^/ur, Uj^Uj/u^ (i = l-3). (3.21) 

With mean field improvement, the one-loop counter-terms of Ur and must be removed from 
perturbative coefficients. 

Here we consider generically the 0{g^) contributions from the mean field to the self-energy. 
From the Feynman rules in Appendix the self-energy and its first derivative with respect to 
P4, on the mass shell (p4,p) = (iMi,0), are 

sWp(^Mi,0) = V^W3rC-5'nWe<'- (3.22) 
sWp(iMi,0) = -g^WAfoa, (3.23) 
SaU^Mi.O) = V4^1e<- (3.24) 



Then, the mean field contribution to the rest mass is 



[1] _ ^QQr [11 , . 

and to the wave function renormalization factor 

which holds for massive and massless (Mg = 0) quarks. The explicit values of u^' and ' 
depend on the definition of the mean field. Since one can easily incorporate the contributions 
from the mean field improvement to the one- loop coefficients, we do not employ a specific 
scheme and quote only the contributions from the loop integrations. 

3.5 Quark bifinear operators 

To obtain improved matrix elements, operators also must be improved [|2^]. As discussed in 
Sec. ^, with the choice r = only the multiplicative factor exp(Mi 07-/2) = \ + ^Miar-\-0{a^) 
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is required at the tree-level. In particular, with r = 1/^ no new dimension- four operator is 
needed to achieve tree-level improvement. At higher loop order the counterpart of the mass- 
dependent factor comes from the quark self-energy through the wave function factor, as seen 



in (3.14), and dimension-four terms are needed. 

Because the tree-level rotation coefficient di vanishes as moflr 
rents of the form 

Jr(x) = ipiTiphix), 



0, we consider here cur- 



(3.27) 



where ipi and iph are the light and heavy quark fields respectively. We consider the vector 
and axial vector currents, so the the 4x4 matrix F is one of 74 (V4), 7j (Vj), 7574 {A4), and 
757j (Aj). We seek the matching factors Zj^ such that Zj^J^ has the same matrix elements 
(for pa„ <C 1) as the continuum bilinear. These matching factors are composed of two parts: 
the wave function of each quark field and the correction to the vertex. Since the former is 
already obtained in previous subsection, here we discuss the vertex corrections. 

The vertex function Ar, which is the sum of one-particle irreducible three-point diagrams, 
can be expanded 

00 

Ar = l+E5^'AM. 

1=1 



(3.28) 



As with the wave function renormalization, Ar is gauge dependent and suffers from infrared 
and mass singularities. For the one-loop term we again subtract of the divergent part. 



A 



[1] 



V^^ + L 



[1] 
r ' 



where, in Feynman gauge, 



L 



[1] 



Cf 
167r2 



(^-^{HG-l)-ln{X^al 



(3.29) 

(3.30) 
-2 for temporal 



The constants are again taken from the continuum expression, so HG 
components of the currents and HG = 2 for spatial components. 

The sought-after matching factor Zj^ is simply the ratio of the lattice and continuum 
radiative corrections: 



■1/2. 7i/2]™>^t 



zli^ArZ 



1/2- 
21 



lat 



(3.31) 



In view of the mass dependence, we write 



il + Moar)'/'Zj, = l + J29'^zl, 



(3.32) 



(=1 



so that the Zj^ are only mildly mass dependent. At the one- loop level we have consistently 
defined the finite lattice parts so that 



(3.33) 



is the desired one- loop coefficient of the matching factor. It is gauge invariant and independent 
of the scheme for regulating the infrared and (light-quark) mass singularities. Numerical results 



for Vp^' and Z^j^ are in Sec. | 
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4 Numerical results of one-loop perturbation theory 



In this section we present our results for the one- loop coefficients. They are obtained numeri- 
cally using the Monte Carlo integration program BASES [ ^ . We give one-loop terms for the 
rest mass, i.e., M^^ and c^^l; for the infrared-finite parts of the wave function renormalization 
factors and the vertex functions, i.e., wf'\ and V^^'; and for the currents' matching fac- 
tors Zj^. In this section we are concerned with zero three- momentum, so for brevity we set 
the temporal lattice spacing a,- = 1- When the spatial lattice spacing is needed, we use the 
anisotropy ^. 

The spatial and the temporal parts of m|^' are listed separately, namely 

g'M^^^^glM^S+glM^}. (4.1) 

so that one could use different (improved) couplings in a practical evaluation of the perturbative 
rest mass. On the other hand, for the other quantities we show the combined values of the 
spatial and the temporal parts, because we are interested mostly in seeing how they behave 
when Mq^ ~ 1 while Mq := MqUt small. 

The spatial and the temporal parts of m|^^ , and m|^^ respectively, are listed in Table |l| 
for a range of Mq < 1 at four values of 1, 2, 3 and 4. We plot J-^^ vs. Mq in Fig. |^ and the 
numerical values are given in Table ^. One sees that the mass dependence is significant, but 
not drastic. 

Table ^ lists the one-loop corrections wj^'^'^ and wj^^ to the massless and massive quark 
wave function renormalization factors. The mass dependence is shown in Fig. ^. Here the 
introduction of anisotropy is seen to reduce the mass dependence greatly. Since wj^^ connects 
smoothly to wj^'^'^ , one sees that we have subtracted the infrared singularities in a consistent way. 



Table 1: The one- loop correction to the rest mass, M[ , for various values of Mq at four values 
of 



Mq 


e = i 


e = 2 


e = 3 


e = 4 


m[]1 0.01 


0.00111(3) 


0.00151(1) 


0.001477(7) 


0.001456(5) 


0.02 


0.00286(3) 


0.00266(1) 


0.002571(7) 


0.002534(5) 


0.05 


0.00624(3) 


0.00543(1) 


0.005116(7) 


0.004958(5) 


0.10 


0.01036(3) 


0.00868(1) 


0.008017(7) 


0.007563(5) 


0.20 


0.01620(3) 


0.01282(1) 


0.011177(7) 


0.009938(5) 


0.30 


0.02002(3) 


0.01504(1) 


0.012449(6) 


0.010507(4) 


0.50 


0.02434(2) 


0.01665(1) 


0.012582(6) 


0.009810(4) 


1.00 


0.02694(2) 


0.015300(9) 


0.009960(4) 


0.006959(3) 


m[^J 0.01 


0.00251(1) 


0.002018(5) 


0.001895(3) 


0.001847(2) 


0.02 


0.00469(1) 


0.003674(6) 


0.003423(3) 


0.003333(3) 


0.05 


0.01049(2) 


0.007864(7) 


0.007236(4) 


0.006959(3) 


0.10 


0.01878(2) 


0.013606(7) 


0.012275(5) 


0.011562(4) 


0.20 


0.03222(2) 


0.022430(8) 


0.019480(5) 


0.017757(4) 


0.30 


0.04305(2) 


0.029095(7) 


0.024581(6) 


0.021753(5) 


0.50 


0.05988(2) 


0.038651(9) 


0.031177(6) 


0.026528(6) 


1.00 


0.08646(2) 


0.05215(1) 


0.039395(7) 


0.031916(6) 
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Table 2: The non- logarithmic part of the mass renormalization factor c^^l. 



Mo 


e = l 


e = 2 


e = 3 


e = 4 


cW 0.01 


0.143(2) 


0.092(1) 


0.0798(4) 


0.0751(3) 


0.02 


0.145(1) 


0.092(1) 


0.0790(2) 


0.0738(2) 


0.05 


0.1438(4) 


0.0900(2) 


0.0760(1) 


0.0690(1) 


0.10 


0.1408(2) 


0.08701(8) 


0.07104(6) 


0.06181(4) 


0.20 


0.1370(1) 


0.08197(5) 


0.06309(3) 


0.05070(2) 


0.30 


0.13374(7) 


0.07828(4) 


0.05744(2) 


0.04354(2) 


0.50 


0.13004(5) 


0.07365(2) 


0.05110(1) 


0.03662(1) 


1.00 


0.12787(3) 


0.07043(1) 


0.04781(1) 


0.03473(1) 



0.2 



0.15 



0.1 



0.05 



- EKM(linear approx. in M^) 

• ^=1 
■ ^=2 

♦ ^=3 
A ^=4 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 3: The non-logarithmic part of the mass renormalization factor c'^J . The dashed line is 
the linear approximation, based on Ref. [ 17|. 
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Table 3: The one-loop correction to the light and the heavy quark wave functions. 



Mn 



0.00 0.08194(2) 0.02994(1) 0.01911(1) 0.01605(1) 



W, 



w 



0.01 
0.02 
0.05 
0.10 
0.20 
0.30 
0.50 
1.00 



0.08009(19) 

0.07887(18) 

0.07537(15) 

0.06949(14) 

0.05892(11) 

0.05072(11) 

0.03833(7) 

0.01908(7) 



0.02913(13) 

0.02892(12) 

0.02774(11) 

0.02649(9) 

0.02417(7) 

0.02258(6) 

0.01979(5) 

0.01388(5) 



0.01890(11) 

0.01918(10) 

0.01929(8) 

0.02009(7) 

0.02137(5) 

0.02209(5) 

0.02233(4) 

0.01715(5) 



0.01621(11) 

0.01634(9) 

0.01765(7) 

0.01981(6) 

0.02354(5) 

0.02561(4) 

0.02594(4) 

0.01856(5) 



0.1 



0.05 



• ^=1 
■ ^=2 

♦ ^=3 

A ^=4 



0.2 0.4 0.6 0.8 



Figure 4: Mass dependence and ^ dependence of the one-loop correction to the heavy quark 
wave function. 
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Table 4: The one-loop correction to the axial vector current vertex function. 





Mn 


f = 1 


f = 2 


f = 3 


£ = 4 




0.00 


0.03449(1) 


0.01005(1) 


0.00033(1) 


-0.00472(1) 




0.01 


0.03425(13) 


0.01028(17) 


0.00071(10) 


-0.00369(10) 




0.02 


0.03399(10) 


0.01061(10) 


0.00134(8) 


-0.00309(8) 




0.05 


0.03426(8) 


0.01145(7) 


0.00337(6) 


-0.00036(6) 




0.10 


0.03435(6) 


0.01294(5) 


0.00628(5) 


0.00399(5) 




0.20 


0.03396(4) 


0.01294(4) 


0.01155(3) 


0.01142(3) 




0.30 


0.03389(4) 


0.01543(4) 


0.01585(3) 


0.01717(3) 




0.50 


0.03337(3) 


0.02160(2) 


0.02230(2) 


0.02479(2) 




1.00 


0.03309(2) 


0.02818(2) 


0.03187(2) 


0.03465(1) 




0.00 


0.03450(1) 


0.02669(1) 


0.02460(1) 


0.02454(1) 


0.01 


0.03428(5) 


0.02629(4) 


0.02391(3) 


0.02359(3) 




0.02 


0.03410(3) 


0.02602(3) 


0.02352(2) 


0.02294(2) 




0.05 


0.03407(2) 


0.02551(2) 


0.02237(2) 


0.02121(2) 




0.10 


0.03378(2) 


0.02463(2) 


0.02075(1) 


0.01896(1) 




0.20 


0.03342(1) 


0.02326(1) 


0.01854(1) 


0.01624(1) 




0.30 


0.03302(1) 


0.02223(1) 


0.01716(1) 


0.01490(1) 




0.50 


0.03253(1) 


0.02093(1) 


0.01582(1) 


0.01397(1) 




1.00 


0.03176(1) 


0.01950(1) 


0.01512(1) 


0.01429(1) 



Tables ^ and ^ list the one- loop corrections to the axial vector and the vector current vertex 
functions, respectively. The mass dependence is shown in Fig. ^ and ^. For the axial vector 
current, the mass dependence with anisotropy is larger than with = 1, but still small. 

These results are combined to obtain the one-loop part of the matching factors according 
to Eq. (Ul). The results are listed in Table ^ and plotted in Figs. |^ and The magnitude 
of the one- loop correction decreases as £ increases, in all cases. Figures ^ and ^ are the most 
important results of this section. They are relevant to phenomenological applications to charm 
physics. Moreover, these results test, at the one-loop level, whether the matching factors are 
well-behaved in the interesting region with small Mq but Mq^ ~ 1. For this reason we have 
plotted them, as in Fig. ||, not against Mq but Mq^. By inspection of Figs. and ^, one can 
see that the small Mq Taylor series continues to be a good approximation to the full mass 
dependence for all Mq^ < 1, for S, > 2. Had we found a stronger mass dependence (like that 
of Ct*^! for r = 1), one would begin to doubt the feasibility of the ideas laid out in Ref. [ Q| also 
for our choice r = l/£. 

For £ = 1, our results should reproduce previous calculations on isotropic lattices. In the 
case of the mass and the wave function renormalization, we (independently) reproduced the 
full mass dependence of Ref. [ |l7|. For the matching factors of the currents, only the result for 
the massless quark is (independently) available [ and we find agreement. 

In conclusion, in the whole region of Mq and £ we surveyed, we found good behavior 
connecting the continuum limit with the region of practical interest. For charmed hadrons, 
the target region of the heavy quark mass Mq is around 0.1-0.3 on lattices with anisotropy 
£=3-4. The required one- loop coefficients of the renormalization factors are easily obtained by 
interpolating the values in the tables using, for example, spline interpolation. 
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Figure 5: Mass dependence and dependence of the one-loop vertex corrections to the axial 
vector current. 
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Figure 6: Mass dependence and ^ dependence of the one-loop vertex corrections to the vector 
current. 
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Figure 7: zj| and zj] vs. Mq^. 
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Figure 8: 



Z^yj and Z^y^ vs. Mqo 
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Table 5: The one-loop correction to the vector current vertex function. 





n nn 
u.uu 


n V\A7A0^^^ \ 


n 09071 ^^\ 


n nnfiQ'^/'i ~\ 


n nnnAQ/ 


1 ^ 
^) 


U.Ui 


U.U4 / 4U(^iO ) 


U.UZUoDl^ii ) 


n nnfi/i 1 n n^ 


n nni 99/ 

— U.UUiZZI 


1 n\ 

iUj 




n no 
u.uz 


U.U4D ( 1 (^iUj 




n nnR/i 1 /q^ 


n nnnQn/ 
— u.uuuyui 






n n^ 




U.UZUUDl^D ) 


n nnfi'^1 


n nnnon/ 
— u.uuuyui 


oj 




n 1 n 




u.uiyo 1 \0) 


n nnfi99/'f^'i 


n nnnfii^/ 

— U.UUUDOI 


7"! 




u.zu 




n ni on/i 
u.uiyu4(^4 j 


n nnfii 7('i\ 


n nnni fi/ 

— U.UUUiOl 






U.oU 




n ni 81 


n nnfii 1 


n nni 1 fi/ 

— U.UUiiOl 






0.50 


0.04165(3) 


0.01669(3) 


0.00579(2) 


-0.00188( 


2) 




1.00 


0.03629(3) 


0.01308(2) 


0.00439(2) 


-0.00179( 


1) 




0.00 


0.04748(1) 


0.04784(1) 


0.04573(1) 


0.04326(1 






0.01 


0.04727(4) 


0.04764(4) 


0.04570(3) 


0.04320(3 






0.02 


0.04746(3) 


0.04797(3) 


0.04589(3) 


0.04366(2 






0.05 


0.04779(2) 


0.04843(2) 


0.04665(2) 


0.04436(2 






0.10 


0.04825(2) 


0.04923(2) 


0.04762(1) 


0.04568(1 






0.20 


0.04917(1) 


0.05075(1) 


0.04967(1) 


0.04807(1 






0.30 


0.05007(1) 


0.05219(1) 


0.05151(1) 


0.05020(1 






0.50 


0.05175(1) 


0.05481(1) 


0.05467(1) 


0.05369(1 






1.00 


0.05510(1) 


0.05977(1) 


0.06046(1) 


0.05980(1 





5 Conclusion 

In this paper, we have studied the 0{a) improvement of Wilson quarks on anisotropic lattices. 
At the tree- level we find that a certain choice of the parameters, r = 1/^ |l^], is well- 
behaved in the region of practical interest for charmed hadrons, namely Moa^ ~ 1, while MqUt 



is small. On the other hand, with a different choice, r = 1 [ ^, 11, 12 1, continuum behavior 
is reached only for Moa^ <C 1. With this latter choice a non-relativistic interpretation [ ^, Q 
is still possible, but a mass-independent renormalization, which was proposed in Ref. [ ^, is 
obstructed. 

The choice r = 1/,^ also simplifies tree-level 0(a) improvement. The action does not re- 
quire separate temporal and spatial hopping parameters. The currents require mass-dependent 
matching factors, but no intrinsically dimension-four terms. 

We therefore have started to examine the behavior of this choice at the one-loop level. We 
have computed the one-loop contributions to the rest mass and to the matching factors of the 
vector and axial vector currents. The matching factors depend significantly on ^. A more 
critical observation is that they are well approximated by Taylor expansions 

Zj^{^,Moar) ^ Zj^iCO) + MoarZ'j^{^,0) (5.1) 

for MqUcj = ^MQttr < 1 and ^ = 2-4. This region encompasses the one suitable for the charmed 
quark with currently available computer resources. 

There are several issues that remain to be studied. The first is to compute the one-loop cor- 
rections to the ratio of hopping parameters C, the clover coefficients cb and c^;, and dimension- 
four terms in the currents. The calculation of C, is especially difficult, because it requires the 
one-loop kinetic mass M^^'. As at the tree-level, it is crucial to compute the full mass depen- 
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Table 6: The one-loop term of the matching factor between the continuum and the lattice 
theories for the axial vector and the vector currents. The spatial and the temporal contributions 
are summed. 



Mo 


e = l 


i = 2 


? = 3 


e = 4 


Z^l[ 0.00 


-0.11643(3) 


-0.04000(3) 


-0.01944(3) 


-0.01133(2) 


0.01 


-0.1153(2) 


-0.0398(2) 


-0.0197(2) 


-0.0124(2) 


0.02 


-0.1144(2) 


-0.0400(2) 


-0.0205(1) 


-0.0131(1) 


0.05 


-0.1129(2) 


-0.0403(1) 


-0.0226(1) 


-0.0165(1) 


0.10 


-0.1101(1) 


-0.0412(1) 


-0.02588(8) 


-0.02192(9) 


0.20 


-0.1044(1) 


-0.04248(8) 


-0.03178(6) 


-0.03122(6) 


0.30 


-0.1002(1) 


-0.04423(8) 


-0.03644(6) 


-0.03800(5) 


0.50 


-0.09351(8) 


-0.04646(6) 


-0.04302(5) 


-0.04579(5) 


1.00 


-0.08360(7) 


-0.05009(5) 


-0.04999(4) 


-0.05196(4) 


Za\ 0-00 


-0.11644(3) 


-0.05664(2) 


-0.04371(2) 


-0.04059(1) 


0.01 


-0.1153(2) 


-0.0558(1) 


-0.04291(9) 


-0.03972(9) 


0.02 


-0.1145(1) 


-0.05545(9) 


-0.04267(8) 


-0.03914(7) 


0.05 


-0.1127(1) 


-0.05436(8) 


-0.04157(6) 


-0.03807(6) 


0.10 


-0.1095(1) 


-0.05285(7) 


-0.04035(5) 


-0.03689(5) 


0.20 


-0.10385(8) 


-0.05032(5) 


-0.03877(4) 


-0.03604(4) 


0.30 


-0.09935(7) 


-0.04849(5) 


-0.03776(4) 


-0.03573(4) 


0.50 


-0.09266(6) 


-0.04580(4) 


-0.03654(4) 


-0.03497(3) 


1.00 


-0.08227(5) 


-0.04142(4) 


-0.03324(4) 


-0.03160(3) 


Z^vl 0.00 


-0.12943(4) 


-0.05066(2) 


-0.02606(2) 


-0.01556(3) 


0.01 


-0.1284(2) 


-0.0499(2) 


-0.0254(2) 


-0.0149(2) 


0.02 


-0.1272(2) 


-0.0496(2) 


-0.0256(1) 


-0.0153(1) 


0.05 


-0.1256(2) 


-0.0489(1) 


-0.0255(1) 


-0.0160(1) 


0.10 


-0.1220(2) 


-0.0478(1) 


-0.02582(8) 


-0.0173(1) 


0.20 


-0.1158(1) 


-0.04610(8) 


-0.02640(7) 


-0.01996(6) 


0.30 


-0.1104(1) 


-0.04442(7) 


-0.02670(6) 


-0.02199(6) 


0.50 


-0.10178(8) 


-0.04156(6) 


-0.02651(5) 


-0.02288(5) 


1.00 


-0.08680(7) 


-0.03499(5) 


-0.02252(5) 


-0.01910(4) 


Z^yl 0.00 


-0.12942(3) 


-0.07779(2) 


-0.06484(2) 


-0.05931(2) 


0.01 


-0.1283(1) 


-0.0772(1) 


-0.06470(9) 


-0.05933(9) 


0.02 


-0.1279(1) 


-0.07740(9) 


-0.06503(8) 


-0.05985(8) 


0.05 


-0.1264(1) 


-0.07728(8) 


-0.06585(7) 


-0.06121(6) 


0.10 


-0.1240(1) 


-0.07744(7) 


-0.06721(5) 


-0.06361(5) 


0.20 


-0.11960(8) 


-0.07781(5) 


-0.06990(4) 


-0.06787(4) 


0.30 


-0.11639(8) 


-0.07846(5) 


-0.07211(4) 


-0.07103(4) 


0.50 


-0.11188(6) 


-0.07968(4) 


-0.07538(4) 


-0.07469(3) 


1.00 


-0.10561(5) 


-0.08168(4) 


-0.07859(4) 


-0.07711(4) 
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dence, so one can check whether low-order Taylor expansions work well for MqUij ~ 1. Only 
with the full mass dependence can one check whether ^, which comes with the couplings in 
the action, and M^ar, which also comes from the on-shell condition, come together to form 
Motta = (,MQar- If not, then one could proceed with a non-perturbative calculation of C(CiO), 
C'(e,0), cb(C,0), ce{^,0), etc. 

A more practical problem is to define renormalized couplings. The scale-setting scheme of 
Brodsky, Lepage, and Mackenzie (BLM) is usually a good way to absorb the dominant part 
of two- and higher-order contributions [ On an anisotropic lattice, it may make sense 

to define separate scales for temporal and spatial gluons. These results are of interest in any 
case: even if anisotropic lattice calculations require a non-relativistic interpretation for heavy 
quarks, anisotropy remains a useful tool for improving the signal-to-noise ratio. 

Finally, after these problems are resolved, it will be important to combine the results with 
numerical simulation data to obtain the matrix elements relevant to experimental measurements 
of charmed hadrons. 
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A Feynman rules 



The Feynman rules for perturbative calculation are shown the same as in Ref. [ 17 1 except for 
two points. One is that the cb is replaced by rgCB as was mentioned in Sec. p.l| and the other is 
the gluon propagator. With a gauge fixing term that is symmetric under exchange of temporal 
and spatial axis, the free propagator of the gauge field, in Feynman gauge, is 

lj = u = 4 

G^,Ak) = { G„{k) ^ = z.<4 (A.l) 




— ^ — — ^ — — ^ — (g) — ^ — 

V T V P a P 

(a) Vr\p) (b) Vr\p) 

Figure 9: Feynman rules required for the one-loop perturbation theory. 
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Gr = ^G^{k). (A.2) 

where we replaced the anisotropy parameter with the tree- level value ^. A fictitious gluon 
mass A is introduced to regulate infrared divergences. 

The Feynman diagram for the counter-term from the mean field is obtained by expanding 

= l + g2^W+0(54), (A.3) 

= l + (A.4) 

With the replacement of the link variables as Eq. ( p.21| ) in the action ( |2.5| ), the Feynman rules 
required for the the one-loop calculation are 

vf^'^\p) = g^i4W[74isinp4 — COSP4] (A. 5) 

V(^^)(p) = 5H''E[C7.-^sinp, -rC cosp,] (A.6) 
i 

with the diagrams (a) and (b) in Fig. 9, respectively. 



B Explicit expressions of one-loop corrections 

In the following, we show the explicit representations of the self-energy and the vertex correc- 
tion. In order to simplify the expressions, we introduce the following abbreviations: 

= cosA;^, = sin/c^, 

Cf, = cos {\ki}j , = sin {\kf^ , 

C4 = cos(iMi + ^4), 5*4 = sin(«Mi + ^4), 

(74 = cos (iMi + iA;4) , 54 = sin (zMi + ^^4) . (B.l) 

3 3 3 

s2 = ^ SiSi, = I] ^^i'^i^ • C] = E sjcl (B.2) 
1=1 i=l 1=1 

To reduce the volume of notation, we also define 

c'b = rCcB, c'e = rcE, (B.3) 

which always appear in these combinations. To reduce the Dirac matrix structure, it is conve- 
nient to introduce = il, with the upper (lower) sign for massive quarks (anti-quarks) on 
the external leg, and G and H, defined by 

GT = 74r74, HT = 7^r7^, (B.4) 

with an implied sum on /x. The following quantities are convenient for representing the one-loop 
expressions below and in our integration programs: 

^9 = -iGc4 - S4, = -iG^Ci - S4, (B.5) 

B'i = -iGsi + 2sl + 2rCs^, = -iG^ Si + 1 + rriQ - Ci + Ir^s^ , (B.6) 

^9 = -iQ + I4GS4, = -iC + yEG^S4, (B.7) 

J? = -iGc4 + S4, = -iG^Ci + Si, (B.8) 
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and 



B'^ = iGs4 + 2sl + 2rCs'^, = iG^ + 1 + mo - + 2rCs'^ , (B.9) 

ji = iGc4 + S4, J'^ =iG^C4 + Si. (B.IO) 

The symbols with superscript q are essentially massless versions of those with superscript Q. 

Prom the Feynman rules in Appendix the contributions to the self-energy from the 
rainbow diagram, Fig. |3|(a), are 



S,(iMi,0) 
EW(iMi,0) = Cf 



(iMi, 0) + (iMi, 0) + 0(/) 



(B.ll) 



(2vr)' 



-SgiiMi + k)G„{k) c?B^{EQf + r'^C'^s^B^ 



+ires^EQ + [iCc'^EQ + 14bQ) (^2^2 _ [s . c])] , (B.12) 



EW(iMi,0) 



(2vr)4 



SQiiMi + k4,k)Gr{k) 



where 

S'q(p) = sin2p4 + C^^siv?Pj + [\{pi + rCp^) +mo} 

j 

Similarly, the contributions from the tadpole diagram. Fig. ^b), are 



<7^Sj,^j(zMi,0) +5^Sj,^J(zMi,0) + 0(5^), 



(B.13) 
(B.14) 

(B.15) 
(B.16) 



i^F J -^0^Gr{k) (G^sinhMi + coshMi) . (B.17) 



The derivative of the self-energy with respective to p4 is separated into temporal and spatial 
contributions 

E(zMi,0) =g2i]W(iMi,0) + 5?sW(iMi,0) + O(/) (B.18) 
The contributions from Fig. 0(a) are 



rriQ 



(B.19) 



j^G„{k) [-25Q(iMi + ki, kfSi (l + 2rCs^ + 

X [c\E'^fBQ + r^^s^B^ + K^s^i^^ + [iCc'^E'^ + ic'|i?«) (s^c^ _ . c])} 
+5q(zMi + k,, k) {c\EQfj<^ + r^C's^jQ + icg(s2g2 _ . -j^Qj 

EW(iMi, 0) = -iC^ J ^^Gr{k) [-2SQ{iNh + k^, kfS^ (l + 2rCs^ + mo) 
X + ^2 (<4c4^«G« + icicl^Q) } 

+5q(zMi + ki, k) [-2i{A'^ fB'iG^ + C4c4s2a« + icgcls^ JQ + J^} 

There is only one contribution from Fig. ^(b) 

d^k 



(B.20) 



sW(iMi,0) = -|C7F 



(2vr)' 



-G^(A;)[GQ cosh Ml + sinhMi], 



(B.21) 
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with a temporal gluon. 

The vertex function is spht as follows: 



(B.22) 



and the contributions are 



A 



[1] 
rcr 



Cf 



SQ{tMl+k4,k)Sq{k)G^{k) 



+ U{H - G) {^(s^c^ - [s ■ c]) + rCs'} {E^B'' + b"^ E"^) 



A 



[1] 



+ K {^4(s2c2 - [s ■ c]) [{H - Gf - 3) + rCs^l (^g^Q + ^g^Q) 
+ {^c'K.^c^ - [. . c])({H - Gf - 3) + r\H'^ B^B^ 

- e {^(s^c^ - \s ■ c\) {{H - Gf -3) + [s. c]} E^E^ 

- ^c'ie^'i^'c' - • mH - G) - lr\'s's\H - G)] T, (B.23) 
Cf 1 -^^SgitMi + ki,k)Sg{k)Gr{k) X 

+ ic4c4s2G) [a^B^ + ICc'eC^s'^G^ 



+^{H - G)s^{iCA'' + ^c'eC^B'g) + ^c'eC^B^G^) 



(B.24) 
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